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O . Abstract 

Abstract telegrapher's equations and some random walks of Poisson type are shown 
to fit into the framework of the Hamiltonian formalism after an appropriate time- 
dependent rescaling of the basic variables has been made. 

7— I . 
> _ 

^ ■ § 1. Introduction 
O 

Time evolution of random processes differs in one essential respect from evolution of 
conservative systems in general and Hamiltonian systems in particular. As a great number 
of various limit theorems attest, the final states loose all the information about the initial 
conditions. Thus, there is nothing to be conserved, no constants of motion can exist, and 



no use can be made of the powerful machine of the Hamiltonian formalism. Or so it seems. 



Sometimes there are things which do not change during time evolution, such as rates of 
decay and other universal exponents. This suggests that one may hope to find constants of 
motion and even Hamiltonian forms in at least some probabilistic systems provided one is 
willing to make time-dependent rescalings of the basic dynamical variables. Besides, there 
is something like a precendent in the history of attempts to quantize dissipative systems, 



^ ■ a close relative of random processes. The simplest of such systems is a particle moving on 

a line under the influence of a harmonic force and a friction: 

x + 2kx + bx = 0, (1.1) 

where: x = x(t), x £ R 1 , is the position of the particle; k and b are constants; and overdot 
denotes the time-derivative. Clearly, for k ^ the equation (1.1) is not a Hamiltonian (or 
a Lagrangian) system as it stands. However, set 

x{t) = X(t) e~ kt . (1.2) 



Then 



x + kx + bx 



X + (b — k ) X e~ kt (1.3) 
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and we get for X the equation 

X + (b - k 2 ) X = (1.4) 
which is a Hamiltonian system with the Hamiltonian 

H(p,X) = P -+ K -^X 2 . 

The main theme of this paper is that some (not all!) random evolutions of Poisson type 
can be put into a Hamiltonian form after an appropriate time-dependent rescaling of the 
basic variables has been made. This rescaling into a Hamiltonian form is, naturally, a 
heuristic principle and not a general theorem. We shall see below how this principle works 
and fails to work for a few representative systems. We start in the next Section with 
one-dimensional random walk. Multidimensional generalizations proceed in two different 
directions: abstract telegraphers's equations, Section 3, or random walks in H d and Z d , 
Section 4. 



§ 2. One-Dimensional Random Walk 

Let us consider a particle which moves on R 1 with a constant speed v, and reverses 
direction according to a Poisson process with intensity a. This model had been proposed 
by G.I. Taylor [1] in an attempt to understand turbulent diffusion. It is convenient to 
work with a discrete situation first, and then pass to the continuous limit; the reader will 
find a lucid analysis in Kac [2] whose treatment I follow. So, suppose we have a lattice 
ZAx. Our particle moves with the speed v in the positive or negative direction; after time 
At = Ax/v it changes direction; with probability 1 — a At the direction stays the same. 
Denote by S n the displacement of the particle after n steps, with the initial step taken in 
the positive direction. Given a function (p = (p(x), define the expectation values 

F+(x):={<p(x + S n )), (2.1a) 

F-(x):={<p(x-S n )). (2.1b) 

Thus, F+(x) (resp. F~(x)) is the expectation value of <p(x) after n steps, when the initial 
direction of the walk starting at x is positive (resp. negative). Considering the n th step 
as once removed from the (n — l) st one, we get in the usual way 

F±(z) = (1 - aA^F^x ± vAt) + aAtF^x ± vAt) (2.2) 
which can be suggestively rewritten as 

Ft(x) - F^jx) _ F^jx ± vAt) - F^jx) 

At At (2.3) 

+a [F^x ± vAt) - F^ix ± vAt)] . 
Passing to the continuous limit, we obtain 

dt 8x 
8F~ 8F~ 
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This is the dynamical system we were after. The point of going through the discrete route 
first is the logical ease of deriving equation (2.2) (and similar equations later on). Now 
comes the rescaling. Set 

F ± = / ± e _at . (2.5) 

Then the system (2.4) becomes 

r Qf+ df+ 



dt dx 



+ af- 



^f_ = _ v ^l + af+ 
dt dx 

and this is patently a Hamiltonian system since it can be written in the form 

d_ ( /+ \ = / vd -a \ ( 8H/5f+ \ 
dt[f- J { a vd ) { SH/Sf- ) 

with 



(2.6) 



(2.7) 



H = ^ 2 - ^ ^ (2.8) 
and with 

d := d/dx. (2.9) 

The matrix 

j vd -o \ 
la vd 



(2.10) 



is skewsymmetric constant-coefficient and is, thus, Hamiltonian. See, e.g., [3], Ch. I, for 
the modern point of view on Hamiltonian formalism; all the Hamiltonian matrices below 
are of this simple kind. Note that the Hamiltonian H (2.8) is the first in the infinite series 

h. =- t^Mll^rc^n, „ £Z+ (2.ii) 

of conserved densities of the system (2.6). Indeed, writing 

hi ~ h 2 (2.12) 

when 

(/ii - h 2 ) e Imd, (2.13) 

we have 

dH n 5H n df+ 5H n df- 
dt ~ 5f+ dt Sf~ dt 

= d 2n (f + ) [vd(f + ) + af~] - d 2n (f~) [-vd(f-) + «/ + ] 

~d 2n (f + )ar-d 2n (f-)af + ~0. 
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Thus, H n is a conserved density. Moreover, it is obvious that all the H^s are in involution; 

{H n ,H N }:=X Hn (H N )~0, Vn,iV eZ+ (2.14) 
where Xh„ is the evolution derivation corresponding to the flow with the Hamiltonian H n : 



f + \ = d_ ( f+ \ = ( vd -a\ 6H n /5f + \ 

f- dt\r ) { a vd { 6H n /5f- ) 

' ' (2.15) 

/ v d 2n + l (f+) + ad 2n {f-) \ 

~ y _ vd 2n+l( f -j +ad 2ny+} J " 

Thus, we have an infinite number of commuting flows with an infinity of commuting 
conserved densities. 

We have considered the simplest possible system. Before moving on to more general 
pastures, it is worthwhile to note that the same equations (2.2) arise for the pair of 
functions, p + (x,t) and p~(x,t), describing the probability of finding the particle at the 
point x at the time t, arriving there from the right (for p + ) or left (for p~) (see [4], Ch. I): 

p ± (x, t + At) = (1 - aAt^ix ± vAt, t) + aAtp T (x ± vt, t). (2.16) 

In this form this equation is easy to generalize for the inhomogeneous case and even for 
the case when the particle is allowed to rest (see [4], Ch. I): 

p ± (x,t + At) 

= [1 - a{x)]p ± {x, t) + [a(x ± vAt) - X(x ± vAt)At]p ± (x ± vAt, t) (2.17) 
+A(x ± vAt)Atp T (x ± vAt, t), 

where A(x) is the local intensity of the Poisson process, and l — a(x) is the local probability 
of resting. Passing to the continuous limit we get 

° ° X (2.18) 

dp — d 

~dt = ~ v ~fa,( a ( x } p ~} + X ( X )(P + ~ P~)- 

If A (formerly a) is not a constant, we cannot renormalize the variables p ± by e _A * since 
t will enter explicitly into the motion equations; the system (2.18) in this case cannot 
be converted into a Hamiltonian form. When, however, A is a constant, even though a 
(formerly 1) is not, a Hamiltonian form is possible. Set 

p ± =p ± e - M . (2.19) 

Then the system (2.18) becomes 

dj)^ d ^ 

= v— (crp + ) + Xp~ 



d ! dx (2.20) 
= -v^-{ap ) + \p + 



dt dx 
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which can be rewritten as 




vd -Xa- 1 \ j 5H/5p + 
Act- 1 vd j \ 5H/5p- 



(2.21) 



with 



h = ^?-(p-)\ (222) 

§ 3. Telegrapher's Equation 

The system (2.4) is 2-component first-order in time. "Now the amazing thing is that these 
two linear equations of first order can be combined into a [single] hyperbolic equation" , 
says Kac ([2], p. 500), and proceeds as follows. Set 

F:=\{F+ + F-), G :=1(F+ - F~) (3.1) 



so that 



OF dG . 

(3.2b) 

whence 

3 2 F OF 2 8 2 F 

^F + 2a m= v ^ (3 - 3) 

which is the telegrapher's equation. Rewritten as 
_ 2 d 2 F 2adF _d 2 F 

v 'W + ^ 2 ~m ~ dx 2 ( ' 

it can be considered as a singular perturbation of the diffusion equation 

1 dF _ d 2 F 
D~dt ~ ~dx 2 

where 



(3.5) 



1 , 2a 
— = hm — 

D v z 

when both a and v tend to infinity. Now, the diffusion equation assumes unlimited 
speeds of microscopic agents, clearly an untenable thesis in view of special relativity. The 
hyperbolic equation (3.4)can be considered then as a sort of relativistic regularization of 
the classical diffusion and heat equations. Let us now look at Hamiltonian properties of 
this equation, but first we generalize it to the form 
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where L is an arbitrary linear selfadjoint operator in arbitrary number of space dimensions: 

L f = L (3.7) 

and e is a constant (considered small if desired). The case 

L = A 2 (3.8) 

where A is a skewadjoint operator: 

A* = -A (3.9) 

is the most direct generalization of differential equations of telegrapher's type to which 
probabilistic interpretation applies [5]; more about this case later on. Set 

u = ue xt (3.10) 

where A is a constant to be specified presently. Since 

e " (4 + 2 °l) eAl - <w + 2(At + o) I + ^ + 2 » A) ' 

choosing 

X = -a/e (3.11) 
we transform equation (3.6) into equation 

e §r = £(II) ' ( 3 - 12 ) 

where 

L:=L+^j (3.13) 

is again a selfadjoint operator. The second-order equation (3.12), written as a first-order 
system 

du „ 

~dt =u 

% (3.14) 
is easily seen to be a canonical Hamiltonian system: 

(3.15) 



d u \ _ 1 W 5H/5u 
di\u)~\-l 0/1 SH/8u 



with 
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(It is in this place that the selfadjointness of L plays a role). Like for the system (2.6), we 
have an infinity of commuting conserved densities for the system (3.14): 



H n = ±uL n (u)-±-uL n+1 (u), 



neZj 



(3.17) 



The alert reader may have noticed that the Hamiltonian form (2.7) of 1-dimensional ran- 
dom walk (2.4) is different from the canonical Hamiltonian form (3.15) of its generalization 
(3.6). How could this have happenend? The ultimate reason is that the system (2.4) is 
more rigid than the scalar second-order equation (3.6): the latter can be written in a mul- 
titude of ways as a 2-component first order system. For example, a direct generalization 
of the Hamiltonian form (2.7) exists for the case when L = A 2 with a skewadjoint A. Then 
equation (3.12) results from the following Hamiltonian system: 



/ ^ \ 

dt 

du 
\~di J 



( 



X(u) + -u 
-X{u) + -u 



X 



at 




(3.18) 



X 



ae 




(3.19) 



where 



X = ±Ae- 1 ' 2 . 



Again, 



l -uX 2n {u)- l -uX 2n {u), n£Z+, 



(3.20) 



(3.21) 



is an infinite commuting set of conserved densities of the system (3.18). In addition, 
the same equation (3.12) results from the following Hamiltonian system, quite different 
from (3.18): 



/ du \ 

dt 
du 

V dt I 



X{u) + -u\ 



e 

X(u) - -u 




(3.22) 




S/5u 
5/5u 



(uu) 



(3.23) 



In this case, an infinity of commuting conserved densities is given by the formula 

H n = uX 2n (u), n G Z+. (3.24) 
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§ 4. Multidimensional Random Walk 

A particle moves in H d with a constant velocity v € {v j}. After each time interval At, 
there is a change of velocity. The change from Vi to Vj has the probability p^, and we 
take 

Pij = Sij +a>ijAt, (4.1) 

with 

J>y=0, Vi (4.2) 

j 

(In the continuous limit, for the set of states {V{} we have a Markov process £(i) with the 
transition probabilities 

Pij(At) = Sij + aijAt + o(At), (4.3) 

but, as in § 1, it is more convenient to start with the discrete picture.) Denoting by 
Fi = Fi(x, t) the probability of finding the particle coming for its velocity change into the 
point X at time t with the velocity Vi, we have, similar to § 1, 

Fi(x,t + At) = Y,PjiFj(.x - ViAt,t). (4.4) 

3 

By virtue of formula (4.1), in the continuous limit we get 

-^ = -{Vi-V){Fi) + Y,ccjiFj (4.5) 

3 

where 

t>i-V:=X>).A (4.6) 

s=l ° Xs 

Set 

F t = fi e- xt . (4.7) 
Then equation (4.5) becomes 

^ = -(tii-V)(/0 + £/fc/i ( 4 -8) 

where 

:= + A5jj (4.9) 
so that the constraint (4.2) turns into 

5>, = A, Vi. (4.10) 
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We are going to analyze the system (4.8), (4.10) from the Hamiltonian point of view. As 
the Hamiltonian we pick 

H = \Y,^? (4.H) 

i 

with some unknown constants Cj's. The constant-coefficient Hamiltonian matrix 

Bij = -Sij-Vi ■ V + Tij (4.12) 

where T = (Tij) is a constant skewsymmetric matrix, reproduces the motion equations 
(4.8) through the Hamiltonian ansatz 

3 
iff 

/3ji = TijCj (no sum on j). (4.14) 

Let us estimate the proportion of Hamiltonian random walks among all of them. The 
dimension of the latter is the dimension of the space of the /3's subject to the condi- 
tions (4.10), but with the understanding that A is at our disposal. Thus, 

Total dim = N 2 - N + 1 (4.15) 

where N is the number of the different /j's (and also the number of the velocities Uj's). 
From (4.10) and (4.14) we get 



flF-SMi) (413) 



A = E^ = E r ^' = (E r 

i i 

so that 



n \ "-3 

i 



' (4-16) 



(or no conditions for A = 0). Thus, we have to look at the dimension of the image of the 
map T x c i — ► (re)* = /?, where 

c := (ci, . . . ,c N y, c := diag (ci, . . . ,c N ), (4.17) 

and where c is a function of T given by formula (4.16). Let us compute this dimension at 
the point in the /3-space which corresponds to 

r »=(-i o) ; (418) 

thus, we assume that N is even: 

N = 2N. (4.19) 
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Then 



and 



c = A 



-1 
1 



1 

1 



(4.20) 



(4.21) 



(This (3 corresponds to a direct sum of one-dimensional random walks.) Let 

r = T + ef (4.22) 

c = c + ec (4.23) 

t = fa + etf (4.24) 

be an infinitesimal change of the objects under consideration (i.e., e 2 = 0). Denote 



r " -* c )• A '- A C '- C - 



Then 

c = A 
and hence 



-1 
1 



Ae diag (71,..., 7iv) 



-Tiv+i 



-7 2 7v 



71 



V 



Tat 



Taking A and C out of A we are left with the matrix 



+ 



-A B 
B f C 



J 



B- diag ^>i'->E B ar) 



which amounts to an arbitrary matrix B subject to the conditions 

E4 = o. Vj. 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 
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Thus, the dimension of the (3's is 

N 2 -N N 2 - N ,— 2 — N ,—2 — s N 2 , 
—2— + + ( N ~N) = 2(N -N) = — -N. (4.31) 

Taking into account our free parameter A, we finally get the dimension of the space of 
Hamiltonian random walks around the point [5q: 

N 2 

Ham dim = — - N + 1 (4.32) 

which is more than half the total dimension (4.15) of the space of random walks. I conclude 
with a few remarks. 

Remark 4.33. Even for d = 1 the random walk model in this section is more general 
than the one-dimensional model considered in § 1, since we allow N{= 2N) different 
velocities rather than two. The case d = 1 is special, having all possible velocities being 
proportional to each other. This fact leads to an existence of an infinity of commuting 
conserved densities for the system (4.8): 

Hn = lY,Ci[r{ft] 2 . (4.34) 

i 

Remark 4.35. In our calculations of dimensions we had no use for the convective terms 
(Vi ■ V)(/j). Had these terms been absent to begin with, e.g., for ^-independent solutions, 
we would have been dealing with a system of ordinary differential equations for the /j's; 
formula (4.32) in this case provides a (very) low bound on the dimension of such systems 
with a particular Hamiltonian representation. Such systems 

F* = F l f3 (4.36) 

where 

F t :=(F 1 ,...,F N ), (3:=(f3ij) 

appear, e.g., for column-sums of the inverse Kolmogorov equation for a Markov process 
with transition probabilities (4.3): 

j t Vij = ^PikPkj (4.37) 

k 

where 

Pkj ■= a k j + XjSkj (4.38) 

and 

Fj:=J2Pij- ( 4 -39) 

i 

Hamiltonian analysis can be applied to the full system (4.37) and also to the direct Kol- 
mogorov equation 

j t Vij = ^PikPk3- (4.40) 
k 
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Remark 4.41. If the set of all possible velocities {Vi} is not discrete but is a continuous 
one, the sum sign in equation (4.8) turns into an integral sign. The Hamiltonian arguments 
undergo a similar minor modification. 

Remark 4.42. If the randomly walking particle has internal degrees of freedom [6], 
equation (4.4) changes into 

F?(x,t + At) = Y,Pj^(x- Vi At,t) (4.43) 

with 

pJT = M*» + ( 4 - 44 ) 

where indices /x, v refer to the internal states. Equation (4.5) then becomes 
with 

J>^ = 0, Vj>. (4.46) 

i,/x 

The Hamiltonian analysis proceeds as before, with the quadratic Hamiltonian now being 

H = \ E ( 4 - 47 ) 

In the simplest one-dimensional case, where = {±v}, and where everything is invari- 
ant with respect to the reflection x t— > — ar, as in § 1, we have 

dF+ OF, 



dt dx 



(4.48) 



where 

(F±)" = 5^ := = := = a^, 

etc. We also have a vector analog of the telegrapher's equation: Set 

F:=±(F + + F_), G:=\(F + -F_) (4.49) 



so that 



rap dG .„ ._ 



(4.50) 



Then 



5? - 25 ^ = w2 £? + (« 2 - 52 + [ «' 5] ) F (4 - 5i) 

where the matrices a and a are subject to the condition 

(l,...,l)(S + a) = (0,...,0). (4.52) 
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